Conservation laws and thermodynamic efficiencies 
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We show that generic systems with a single relevant conserved quantity reach the Carnot efficiency 
in the thermodynamic limit. Such a general result is illustrated by means of a diatomic chain of 
hard-point elastically colliding particles where the total momentum is the only relevant conserved 
quantity. 
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Conservation laws strongly affect transport properties. 
Conserved quantities may lead to time correlations not 
decaying with time, so that transport is not diffusive 
and is described, within the linear response theory, by 
diverging transport coefficients. This ideal conducting 
(ballistic) behavior can be firmly established as a conse- 
quence of an inequality by Mazur [H-Q which, for a sys- 
tem of size A characterized by M conserved quantities 
Q n , n = 1,---,M, bounds the time-averaged current- 
current correlation functions as 
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where (• • -)t denotes the thermodynamic average at tem- 
perature T. The constants of motion, Q n , are orthogonal 
to each other, i.e., (Q n QmjT = (Qn)TO~n,m, and relevant, 
that is, (JQu)t 7^ for all n. A non-zero right-hand 
side in Eq. ([T]) at the thermodynamic limit implies a fi- 
nite Drude weight for the current J, which in turn indi- 
cates ballistic transport 4, 5]. The impact of motion con- 
stants on the electric and thermal conductivities has been 
widely investigated 0^3|- In particular, anomalous heat 
transport has been discussed for momentum conserving 
interacting systems in low dimensions [1, However, to 
the best of our knowledge, conservation laws have never 
been discussed for coupled flows, in particular in relation 
to the problem of optimizing thermodynamic efficiencies. 

The search of a new technology capable of reducing the 
environmental impact of electrical power generation and 
refrigeration has aroused great interest in thermoelectric- 
ity, namely the possibility to build a type of solid-state 
heat engine capable of converting heat into electricity, or 
alternatively electricity into cooling I3"l3l- The main 
difficulty is connected to the low efficiency of such heat 
engine. We recall that the maximum thermoelectric effi- 
ciency as well as the efficiency at maximum power (l5r - 
[l9j are determined, within the linear response regime and 
for systems with time- reversal symmetry [loj, by the so 
called figure of merit ZT, which is a dimensionless quan- 
tity, a combination of the three main transport properties 



of a material: the thermal conductivity k, the electrical 
conductivity a and the thermopower (Seebeck coefficient) 
S, as well as of the absolute temperature T: 
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T. 



The maximum efficiency is given by 
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where rjc is the Carnot efficiency, while the efficiency 
^(W'max) at maximum output power reads (l5j 
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The only restriction imposed by thermodynamics is 
ZT > 0, so that both efficiencies are monotonous grow- 
ing functions of the figure of merit and 7/ max — » rjc, 
W(W max ) —> when ZT — > oo. It has been suggested 
that the value ZT = 3 is the target to be achieved in 
order to make thermoelectric engines economically com- 
petitive. In spite of recent progress in material science, 
present technology is limited to low ZT materials and no 
clear path has been identified in order to increase effi- 
ciency. 

A promising new approach, based on the theory of dy- 



namical systems, has been recently introduced [22, 
The hope is that the analysis of idealized models may lead 
to some insight on the microscopic mechanisms which 
lead to high figure of merit in more realistic materials. 
While for non-interacting systems, even in the classical 
framework, the energy filtering mechanism [3-0] has 
been shown to allow to reach Carnot efficiency, very little 
is known for interacting particles. Recent numerical and 
empirical evidence has shown that for a one-dimensional 
diatomic disordered chain of hard-point elastically collid- 
ing particles, the figure of merit ZT diverges as the num- 
ber of particles increases [53, 0] . Since it has been ver- 
ified that the energy filtering mechanism does not work 



2 



here [29( , it follows that the divergence of ZT in the ther- 
modynamic limit rests on a different, unknown property. 

In the present paper we analyze and solve this problem. 
Indeed, we show that for systems having a single relevant 
constant of motion, the electric conductivity is ballistic, 
i.e., a oc A, the heat conductivity is subballistic, k oc A" 
with a < 1, and the thermopower is size independent, 
S oc A , so that the figure of merit ZT oc A 1- " oo 
in the thermodynamic limit A — > oo. Our findings are 
illustrated by the above mentioned prototype model of 
interacting one-dimensional system: a diatomic chain of 
hard-point elastically colliding particles, where the total 
momentum is the only relevant constant of motion. 

We start from the equations connecting fluxes and 
thermodynamic forces within linear irreversible thermo- 
dynamics 30l. 31 1 : 



Jp — LppXi + L pu X2, 
Ju = L up X\ + L uu Xi 1 
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where J p and J, 



are the particle and energy currents, 
and the thermodynamic forces X\ = — V(/3/i), X 2 = V/3, 
with p chemical potential and j3 = 1/T inverse temper- 
ature. (We set the Boltzmann constant ks = 1.) The 
Onsager coefficients (i,j = p,u) are related to the 
familiar transport coefficients as follows: 
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where L denotes the Onsager matrix with matrix ele- 
ments Lij and we have set the electric charge of each 
particle e — 1. Thermodynamics imposes det_L > 0, 
Lp P > 0, L uu > 0, and L up = L pu . The figure of merit 



J pp 
reads 



ZT 



i.L U p ~ pL p 
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It diverges (thus leading to maximum efficiency) iff the 
Onsager matrix L is ill-conditioned, that is, in the so- 
called strong-coupling condition, for which the energy 
and particle currents are proportional, J u oc J p , the pro- 
portionality factor being independent of the values of the 
applied thermodynamic forces. 

The Green-Kubo formula expresses the Onsager coeffi- 
cients in terms of correlation functions of the correspond- 
ing current operators, calculated at thermodynamic equi- 
librium 32j, |33j : 



where 
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The real part of can be decomposed into a 5- 

function at zero frequency defining the generalized Drude 
weight T>ij (for i = j = p this is the conventional Drude 
weight) and a regular part L™ s (cj): 



Reiy(w) = 2nV ij S(uj) + Lj?(w). 
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Non-zero Drude weights, ^ for i,j = p,u are 



a signature of ballistic transport, namely Li 



A at 



the thermodynamic limit, and therefore the thermopower 
S = Lp U /(TL pp ) — p/T oc A . 

We now discuss the influence of conserved quantities 
on the figure of merit ZT. We make use of Suzuki's 
formula [2|] for the currents J p and J Ul which generalizes 
Mazur's inequality {1} by stating that, for a system of 
finite size A, 
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where the summation is extended over all the M orthog- 
onal constants of motion which are relevant for the con- 
sidered flows, that is, non-orthogonal to the currents J p 
and J u , i.e., (J p Q n )T ^ and {J u Qn)T ^ 0. (Irrelevant 
constants of motion are not included in the summation 
since they would give zero contribution.) The presence 
of relevant conservation laws implies that the finite-size 
generalized Drude weights 



L>ij (A) 
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are different from zero. If at the thermodynamic limit 
the generalized Drude weight 



V, 



lim lim — — 
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di'(J l (t , )J,(0))T (13) 



is non-zero, then we can conclude that transport is bal- 
listic 34]. Note that in Eq. (fl"3|) the thermodynamic 
limit A — > oo must be taken before the long-time limit 
t — > oo. The below developed theory only applies to 
the cases in which the two limits commute, that is, 
2?tf=:iimA-K»A*(A) 0. 

If there is a single relevant constant of motion, M = 1, 
due to Suzuki's formula (fTTj) (and assuming that the two 
limits A — > oo and t — > commute), the ballistic con- 
tribution to det L vanishes, since it is proportional to 
DppVuu — L) 2 pu , which is zero from (fTTj) . Hence, detL 
grows due to the contributions involving the regular part 
in Eq. (fTO)) , i.e., slower than A 2 , thus implying that the 
thermal conductivity k oc det L/L pp grows subballisti- 
cally. That is, n oc A Q , with a < 1. Since a oc L pp is 



ballistic and S — L pu /TLp P - 
that ZT = <jS 2 T/k oc A 1 "" 



-p/T oc A , we can conclude 
— >■ oo when A — > oo. 
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The situation is drastically different if M > 1 . In this 
case, due to the Schwartz inequality, 



D D — T) 1 = Ha; I 
where 



(a; p ,a; M )>0, (14) 



l/T 



M/T 



(15) 

= J2k=i x P kXuk - The ec l ua lity arises only 
in the exceptional case when the vectors x p and x u are 
parallel. Hence, for M > 1 we expect, in general, det Z cx 
A 2 , so that heat transport is ballistic and ZT cx A . 

In order to illustrate the above general ideas, we con- 
sider a one-dimensional, diatomic disordered chain of N 
hard-point elastically colliding particles with randomly 
distributed coordinates Zi € [0,A], velocities Vi, and 
masses € -{Vi,^} • The numerically observed di- 
vergence of the figure of merit ZT for v\ 7^ vn (27l - E9j | 
can be understood now in terms of the above developed 
theory. Indeed, in this system there is a single relevant 
constant of motion Q± = P, where P = YliLi m i v i is the 
overall momentum [361 ] . In this case the particle current 
J = T^Li Vi and the energy current J u = X}j=i \ m i v f 



Note that the mass current, J m = J2f=i m i v i' equals the 
total momentum P and therefore does not decay, while on 
the other hand J m ~ rhJ p , where m is the average mass 
per particle, hence we expect that (J p (i)J p (0)) does not 
decay either, so that a ~ A. 

It is easy to compute analytically the time-averaged 
correlation functions Cy(A) (from the second line of 
Eq. (fTTjl ) and then the finite-size generalized Drude 
weights 
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Here iV = N% + N2 with N± and the number of parti- 
cles with mass v\ and respectively. Note that, as ex- 
pected from the above theory, D pp (A)D uu (A) — D pu (A) = 
for any system size A. 

To numerically confirm the above results, we com- 
pute the autocorrelation functions of J p and J u , and 
the cross correlation function between them. In doing 
so we apply periodic boundary conditions and assign to 
Ni (N2) particles of mass 1^1 (y 2 ) random initial po- 
sitions and random initial velocities derived from the 
Maxwcll-Boltzmann distribution corresponding to tem- 
perature T. We then evolve the system and compute 
Cij{A,t) = f f*dt'{Ji(t')Jj(Q)) T (for i,j = p,u) up 
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FIG. 1: (Color online) The current-current correlation func- 
tions for various system sizes with A = 256 (red dashed 
curve), 512 (blue dash-dotted curve), and 1024 (black solid 
curve). The temperature T = 1, particle masses are v\ = 1 
and ^2 = v5- It is seen that all the current-current correlation 
functions approach a finite nonzero value as the correlation 
time increases and that the characteristic time scale to ap- 
proach such value is independent of the system size. Note that 
in all the figures the particle density is fixed to be N/A — 1 
and N!=N 2 = N/2. 



to a time t sufficiently long to obtain stable averages, 
thus estimating dj(A) = limt->oo c%j (A, t). We show 
in Fig. Q] that the current-current correlation functions 
(Ji(t)Jj(0))T do not decay to zero as the correlation 
time is increased, implying that Cij(A,t) do not decay 
either and thus indicating ballistic transport. We finally 
estimate the finite-size generalized Drude weights from 
the time-averaged correlation functions as (A, t)/(2A), 
with a sufficiently long time t to approximate the asymp- 
totic value Dij(A) = Cy(A)/(2A). As shown in Fig. O 
the numerically determined values are in very good agree- 
ment with the theoretical values Dij given by Eq. (fT6| . 

Fig. [T] provides clear evidence that the convergence of 
the correlation functions Cij(A, i) to their asymptotic val- 
ues Cij (A) takes place on a time scale independent of the 
system size A. Therefore, Fig. [T] provides a strong indi- 
cation that for the model under investigation the two 
limits A — > 00 and t — > 00 do commute, so that we can 
compute the generalized Drude weights at the thermo- 
dynamic limit as T>ij = lim.A->.oo Dij(A). Note that in 
taking the thermodynamic limit, we keep constant the 
particle density N/A and the ratio Ni/N 2 - 

Finally, we perform a nonequilibrium calculation of the 
various transport coefficients. (Fo r technical details of 
numerical simulations see Ref. 28(.) According to our 
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FIG. 2: (Color online) Comparison between the numerically 
determined finite-size generalized Drude weights (symbols) 
and the analytical values for Dij (lines) given in Eq. (|16[) 
for A = 256. We set T = 1 and z/i = 1 in (a) and v\ = 1, 
^2 = v 7 ^ in (b). 



theory, we expect that all the Onsager coefficients grow 
linearly with the system size in the thermodynamic limit, 
since the generalized Drude weights T>ij are all different 
from zero. This expectation is confirmed by the data 
shown in Fig. [3] Finally, in Fig. |4]we show the transport 
coefficients a, S, k, and the thermoelectric figure of merit 
ZT as a function of the system size. In agreement with 
our theory, we observe that a cx A, while S saturates to 
the value predicted from theory for ballistic transport, 
S = ± -/i) = | , « oc A" with a « 1/3 [13, H|, 
and the growth of the figure of merit in good agreement 
with the dependence ZT oc A 1- ". 

Note that the above conclusions for the thermal con- 
ductivity k and the figure of merit ZT do not hold in the 
integrable case v\ = i>2, where M > 1 since all moments 
of the momentum distribution are conserved quantities. 
In this case k oc A and it is easy to analytically compute 
ZT = 1 Q. 

To summarize, we have shown that for systems with a 
single relevant constant of motion, the thermoelectric fig- 
ure of merit diverges as the system size increases, so that 
the Carnot efficiency is achieved at the thermodynamic 
limit. Such a result has been illustrated in the case of a 
chain of hard-point elastically colliding particles, with a 
remarkable agreement between analytical results, equilib- 
rium and out-of-equilibrium numerical simulations. We 
would like to point out that, while our illustrative model 
is one-dimensional, there are no dimensionality restric- 
tions in our theory, so that it should apply also to two- 
and three-dimensional systems in which total momentum 
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FIG. 3: (Color online) Dependence of the Onsager coefficients 
Lij on the system size A. The straight lines are the best linear 
fitting, Lij oc A. In these nonequilibrium simulations we set 
the temperature T = 1, the chemical potential /1 = 0, and the 
masses v\ = 1 and V2 = \/T9. 
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FIG. 4: (Color online) Dependence of the transport coeffi- 
cients a, S, « and ZT on the system size A, for the same 
parameter values as in Fig. [3] Dashed lines are drawn for 
reference. 



is the only relevant constant of motion. Therefore, our 
paper unveils a rather generic mechanism for increasing 
thermoelectric efficiency in interacting systems. 
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